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1.0  BACKGROUND  INFORMATION 


1.1  Introduction 

When  radioactive  material  is  released  to  the  atmosphere,  estimates  of  radiation  doses 
received  by  persons  downwind  are  generally  made  by  use  of  standard  atmospheric 
dispersion  equations  (ref.  I,  pp.  97-116).  This  procedure  is  adequate  as  a  means  of 
estimating  inhalation  doses,  but  it  may  result  in  very  large  errors  when  evaluating 
external  gamma  doses  from  radioactive  clouds  (2).  Therefore,  a  number  of  cloud 
gamma  dose  models  have  been  developed  (3,  4). 

Most  published  cloud  gamma  dose  models  were  developed  as  part  of  computer  codes 
used  to  estimate  doses  resulting  from  releases  by  nuclear  power  reactors.  Such 
models  may,  when  applied  to  small  radiological  research  facilities,  yield  estimates 
with  large  and  unpredictable  errors.  The  model  described  in  this  study  was  developed 
to  provide  dose  estimates  for  such  releases.  This  model,  which  uses  exposure  rate 
source  geometries  consisting  of  coneentrie  sets  of  hollow  cylinders,  has  been 
designated  the  concentric  cylinder  set  (CCS)  model.  A  more  detailed  discussion  of 
the  background  of  the  CCS  model  has  been  presented  elsewhere  (ref.  5,  pp.  2-7, 
13-19). 

1.2  Purpose 

The  CCS  eloud  gamma  dose  model  was  developed  solely  for  the  purpose  of  providing 
external  gamma  dose  estimates  for  a  variety  of  radioactive  cloud  geometries.  The 
CCS  model  may  be  applied  to  the  simple  plume  geometries  and  long-term  "average" 
meteorological  conditions  postulated  in  most  models.  However,  the  primary  purpose 
of  developing  the  CCS  model  is  to  estimate  doses  from  irregular  clouds. 

An  irregular  cloud  may  possess  any  of  the  following  characteristics:  (a)  finite  size, 
(b)  more  than  one  plume  centerline,  (c)  more  than  one  plume  orientation,  (d)  not 
accurately  described  with  single-valued  "effective"  meteorological  parameters, 
(e)  widely  varying  concentrations  of  radioactivity  within  one  gamma  ray  mean  free 
path  of  the  receptor,  (f)  plumes  resulting  from  widely  vurying  release  rates,  and 
(g)  receptors  at  downwind  distances  where  entrainment  of  released  gases  has  not  yet 
been  completed. 

1.3  Model  Comparisons 

1'he  CCS  model  has  been  compared  to  a  number  of  published  cloud  gamma  dose 
models  for  those  conditions  to  which  the  published  models  can  bo  applied  (ref. 
pp.  78-83).  Comparison  was  also  made  between  calculated  (CCS  model)  exposure 
rates  and  those  obtained  by  measurements  taken  during  the  same  stack  releases  (ref. 
.I,  pp.  113-117).  The  CCS  model  generally  agreed  with  published  models  to  within 
3()"o  when  identical  input  parameters  could  be  used.  The  model  generally  agreed  with 
field  measurements  within  a  factor  of  2. 

A  point-by-point  model  was  developer)  solely  as  a  comparison  model.  The  point-by¬ 
point  model  is  basically  a  computer  code  derived  from  theoretical  eloud  gamma  dose 
equations,  without  the  special  assumptions  used  in  developing  the  CCS  model  (6). 
Only  a  lew  sets  of  release  conditions  were  evaluated  in  this  comparison,  due  to  the 
large  amount  of  machine  time  required  for  this  model.  Comparison  to  the  CCS  model 
showed  consistent  agreement  within  2U'\>. 
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2.0  BASIC  MODEL 


2.1  General  Description 

The  concentric  cylinder  set  model  requires  little  expenditure  of  money  or  man¬ 
hours.  It  is  not  a  computer  code  and  requires  no  computer  for  its  application.  The 
model  is  flexible,  due  to  its  modular  data  sets  and  absence  of  any  decision  making  in 
the  model  proper.  All  information  needed  to  apply  the  model  is  provided  in  this 
report,  except  for  specific  radioactive  source  data.  Values  for  source  gamma  ener¬ 
gies  and  half-lives  may  be  obtained  from  standard  references  (7). 

A  number  of  assumptions  have  been  made  in  developing  the  CCS  model,  particularly 
in  its  application  to  irregular  clouds.  These  assumptions,  which  have  been  discussed 
elsewhere  (ref.  5,  pp.  57,  58,  92),  may  be  summarized  as  follows: 

(a)  Any  radioactive  cloud  may  be  considered,  for  purposes  of  external 
radiation  dose  evaluation,  to  consist  of  a  few  (generally  less  than  20) 
smaller  source  volumes.  Each  source  volume  may  be  described  by  a  set  of 
single-valued  meteorological  and  geometric  parameters. 

(b)  A  Gaussian  diffusion  model  (8)  may  be  used  to  approximate  the  inter¬ 
nal  source  distribution  of  each  source  volume,  or  plume  segment.  This 
condition  holds,  even  if  the  segment  has  undergone  changes  in  direction 
during  the  period  of  its  travel  from  release  point  to  exposure  of  a  recep¬ 
tor.  This  assumption  is  acceptable  for  external  gamma  dose  estimates  but 
not  for  internal  dose  estimates. 

(c)  Assumptions  generally  used  in  cloud  gamma  dose  models  (ref.  1,  pp. 
337-355)  (9)  may  be  used  for  the  CCS  model  when  used  to  evaluate  both 
simple  plumes  and  irregular  clouds,  with  proper  adjustments. 

2.2  Simple  Flume  Geometry 

The  most  useful  application  for  the  CCS  model  is  the  evaluation  of  doses  from  irreg¬ 
ular  clouds.  However,  application  to  simple  plumes  provides  a  better  background  for 
learning  to  apply  the  model,  and  will  be  discussed  first.  Figure  1A  is  a  side  view  of 
an  idealized  simple  elevated  plume,  as  represented  in  the  CCS  model.  The  plume  is 
represented  by  a  train  of  end-to-end  concentric  sets  of  hollow  cylinders.  A  per¬ 
spective  of  several  cylinder  sets  is  shown  in  Figure  IB.  The  plume  centerline  is  taken 
as  the  +  X  coordinate  in  standard  dispersion  calculations. 

Each  cylinder  set  has  a  Gaussian  source  distribution.  The  maximum  cylinder  radius 
included  in  the  dose  summations  is  denoted  sInax,  and  is  equal  to  3  ',  where  o  is 
magnitude  of  the  dispersion  coefficient  (ref.  1.  pp.  100-105).  Although  this  param¬ 
eter  increases  monotonically  and  continuously  in  standard  dispersion  calculations,  it 
is  treated  as  a  step-function  in  the  CCS  model. 

The  minimum  distance  between  the  receptor  and  the  plume  centerline  is  designated 
the  centerline-receptor  distance,  denoted  by  a.  The  length  of  a  set  of  cylinders  is 
denoted  bv  L.  If  a  receptor  is  located  directly  below  the  plume  centerline,  then  a  = 
li,  the  effective  stack  height.  Since  the  CCS  model  represents  a  theoretically 
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conical  geometry  bv  a  set  of  end-to-end  cylinders,  the  two  geometries  become  ident¬ 
ical  only  as  the  cylinder  lengths  approach  zero.  In  practice,  however,  the  gamma 
doses  assoc ia toil  with  the  two  geometries  are  nearly  equal  for  the  conditions  used  in 
this  report. 


2.3  Basic  Concentric  Cylinder  Set  Equation 

The  basic  CCS  exposure  rate  module  is  based  on  the  use  of  a  hollow  cylinder  as  the 
smallest  separate  source  element.  An  integral  solution  exists  for  the  fluence  rate 
from  such  a  source  geometry  (ref.  10,  pp.  397-415).  If  attenuation  is  assumed  to  have 
a  lesser  effect  on  the  exposure,  the  solution  can  be  used  to  approximate  the  expected 
exposure  rate.  In  the  equation  below  used  to  develop  the  exposure  rate  curves, 
integration  is  applied  over  the  cylinder  internal  angle  and  summation  over  the  other 
dimensions  (cylinder  length  and  radii). 

The  exposure  rate  module  equation  for  exposure  rate  is: 


E.R.  = 


TQ'FCA 


(1  +  k r-j)  exp  (-Atrjj  -  s?/2  a2) 

f7~2  2  dh2  T~2  2 

v(Sj  +  n  +  Z.)  -  4  Sj  a 


AH,  with: 


=  exposure  rate  at  receptor,  mR/h 

-  length  of  each  cylinder  in  the  set,  m 
=  radioactivity  release  rate,  mCi/s 

=  radius  of  the  j*-11  cylinder  in  the  set,  m 
=  point  source  gamma  exposure  rate  constant,  R-m^/Ci  -  h 
=  distance  along  the  cylinder  axis,  to  i^  increment,  m 
=  average,  or  effective,  wind  speed,  m/s 

-  product  of  dispersion  coefficients,  m^,  witii  =  ,/ciy  oz 
=  gamma  buildup  constant  in  air,  no  units 

=  gamma  attenuation  coefficient  in  air,  m~* 

-  centerline-receptor  distance,  in 

=  radioactive  decay  correction  factor,  no  units 
=  cylinder  eccentricity  correction  factor,  no  units 
=  air-ground  scattered  gamma  correction  factor,  no  units 

-  distance  between  receptor  and  cloud  volume  corresponding  to  f.p  sj,  m 
=  increment  of  cylinder  length,  along  axis,  m,  and 

=  increment  of  cylinder  radius  in  a  set,  m. 
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2.4  Dispersion  Coeffieient 

Meteorological  parameters  can  be  obtained  by  use  of  a  wind  speed-  and  direction¬ 
recording  system  and  by  visual  observation  of  the  sky.  Values  of  the  dispersion 
coefficient,  o,  are  provided  in  Figure  2.  They  should  be  applied  with  care  in  order 
not  to  be  used  outside  the  limiting  conditions  of  their  applicability  (11). 


Stability  classes  (or  types)  chosen  for  this  study  are  listed  in  the  legend  for  Figure  2. 
They  are  recommended  by  the  U.S.  Nuclear  Regulatory  Commission  (USNRC).  Type 
D  should  be  used  for  all  heavily  overcast  conditions.  The  appropriate  classes  can  be 
ascertained  by  referring  to  Table  1,  modified  as  necessary  from  standard  references 
(12). 


Table  I .  Categories*  of  Meteorological  Stability 


Daytime  Brightness 

Nighttime 

Cloudiness 

Surface  Wind 
Speed  (m/s) 

Strong 

Moderate 

Slight 

50%  or  More 

40%  or  Less 

<2 

A 

A-B 

11 

(F)  + 

(F-G) 

2 

A-B 

B 

C 

E 

F 

4 

B 

B-C 

c 

D 

E 

6 

c 

C-D 

u 

U 

u 

>6 

C 

D 

D 

D 

D 

*  Categories  taken  from  references  11  and  12 

+  Parentfieses  indicate  that  tnese  were  not  included  in  original  typing  scheme. 


3.0  APPLICATION 
3 . 1  Cylinder  Eccentricity 

In  the  basic  CCS  equation,  a  circular  cylinder  geometry  is  assumed;  this  corresponds 
to  isotropic  diffusion  in  the  yz-plane.  This  assumption,  while  not  consistent  with 
observed  diffusion  patterns,  generally  results  in  acceptably  small  errors  when  used 
only  for  external  gamma  dose  estimates.  However,  there  are  cases  in  which  the 
resulting  error  may  be  significant.  So  it  is  necessary  to  include  a  correction  factor, 
C,  for  cylinder  eccentricity. 

Figure  3A  qualitatively  illustrates  the  effect  of  cylinder  eccentricity  on  source- 
receptor  geometry.  Table  2  lists  values  of  the  ratios  of  crosswind-to-vertical  disper¬ 
sion  coefficients  for  the  range  of  conditions  for  which  the  COS  model  might  be 
needed.  Table  3  lists  correction  factors.  Correction  factors  for  cylinder  eccentricity 
exposure  rate  were  obtained  by  FORTRAN  summation  over  cylinder  set  volumes. 
Values  of  C  arc  listed  for  the  following  parameter  values:  ratio  of  iy/oz.  ratio  of 
-/a,  ratio  of  L/a.  and  values  of  the  angle  which  is  the  angle  between  the  line  a 
and  the  vertical,  l  or  a  ground  level  receptor  and  with  ov  >  oz.  C  is  always  less 
than  1.0,  and  can  be  ignored  for  conservatism. 
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Table  3.  Correction  Factors  for  Cylinder  Eccentricity 


4 

(o) 

a/a 

L/a 

Factor, 

C,  for 

0  /o 

y  z 

= 

4 

(o> 

o/a 

L/a 

Factor, 

C,  for 

yv 

0.5 

1.0 

2.5 

4.0 

6.5 

0.5 

1.0 

2.5 

4.0 

6.5 

0 

0.25 

<  1 

1.4 

1 

0.8 

0.6 

0.5 

45 

0.25 

<  1 

1.0 

1 

0.9 

0.8 

0.7 

1-2 

1.2 

1 

0.8 

0.7 

0.6 

1-2 

1.0 

1 

0.9 

0.8 

0.7 

>  2 

1.2 

1 

0.8 

0.7 

0.6 

>  2 

1.0 

1 

0.9 

0.8 

0.7 

0.75 

<  1 

1.2 

1 

0.6 

0.4 

0.3 

0.75 

<  1 

0.8 

1 

0.9 

0.6 

0.4 

1-2 

1.1 

1 

0.7 

0.5 

0.4 

1-2 

0.8 

1 

0.9 

0.7 

0.5 

>  2 

1.1 

1 

0.7 

0.6 

0.4 

>  2 

0.8 

1 

0.9 

0.7 

0.5 

1.25 

<  l 

1.0 

1 

0.8 

0.6 

0.5 

1.25 

s  I 

0.9 

1 

0.9 

0.8 

0.6 

1-2 

1.0 

1 

0.9 

0.7 

0.6 

1-2 

0.9 

1 

1.0 

0.8 

0.7 

>  2 

1.0 

1 

0.9 

0.7 

0.6 

>  2 

0.9 

1 

1.0 

0.9 

0.7 

1.75 

<  1 

0.9 

l 

0.9 

0.8 

0.6 

1.75 

<  1 

0.9 

1 

1.0 

0.9 

0.7 

1-2 

0.9 

1 

0.9 

0.8 

0.7 

1-2 

0.9 

1 

1.0 

0.9 

0.8 

>  2 

0.9 

1 

0.9 

0.9 

0.7 

>  2 

0.9 

1 

1.0 

0.9 

0.8 

2.25 

1 

0.9 

1 

0.9 

0.8 

0.7 

2.25 

<  1 

0.9 

1 

1.0 

0.9 

0.8 

1-2 

0.9 

1 

1.0 

0.9 

0.8 

1-2 

0.9 

1 

1.0 

1.0 

0.9 

-  2 

0.9 

1 

1.0 

0.9 

0.8 

-  2 

0.9 

1 

1.0 

1.0 

0.9 

2.75 

s  1 

0.9 

1 

1.0 

0.9 

0.8 

2.75 

<  1 

0.9 

1 

1.0 

1.0 

0.9 

i  -2 

0.9 

1 

1.0 

1.0 

0.9 

1-2 

0.9 

1 

1.0 

1.0 

0.9 

>  2 

0.9 

1 

1.0 

1.0 

0.9 

>  2 

0.9 

1 

1.0 

1.0 

0.9 

30 

0.25 

1 

1.1 

1 

0.8 

0.7 

0.6 

90 

0.25 

v  1 

0.8 

1 

1.1 

1.0 

0.9 

1-2 

1.1 

1 

0.9 

0.7 

0.6 

1-2 

0.9 

1 

1.1 

1.0 

0.9 

■  2 

1.1 

1 

0.9 

0.8 

0.6 

>2 

0.9 

1 

1 .0 

1.0 

0.9 

0.75 

<  1 

1.0 

1 

0.6 

0.4 

0.3 

0.75 

<  1 

0.6 

1 

1.3 

1.4 

1.4 

1-2 

1.0 

1 

0.7 

0.5 

0.4 

1-2 

0.7 

1 

1.3 

1.2 

1.2 

2 

1.0 

1 

0.7 

0.5 

0.4 

>  2 

0.7 

1 

1.2 

1.2 

1.2 

1.25 

<  1 

1.0 

1 

0.9 

0.7 

0.5 

1.25 

<  1 

0.8 

1 

1.2 

1.1 

1.1 

1-2 

1.0 

1 

0.9 

0.7 

0.5 

1-2 

0.8 

1 

1.1 

1.1 

1.1 

>  2 

1.0 

1 

0.9 

0.7 

0.6 

>  2 

0.8 

1 

1.1 

1.1 

1.0 

1.75 

<-  1 

0.9 

1 

0.9 

0.8 

0.6 

1.75 

<  1 

0.8 

1 

1.1 

1.1 

1.1 

1-2 

0.9 

1 

0.9 

0.8 

0.7 

1-2 

0.8 

1 

1.1 

1.1 

1.0 

>  2 

0.9 

1 

1.0 

0.9 

0.7 

'  2 

0.8 

1 

1.1 

1 . 1 

1.0 

2.25 

<  1 

0.9 

1 

1.0 

0.9 

0.8 

2.25 

-  1 

0.9 

1 

1.1 

1  .1 

1.0 

1-2 

0.9 

1 

1.0 

0.9 

0.8 

1-2 

0.9 

1 

1.1 

1.1 

1.0 

>  2 

0.9 

1 

1.0 

0.9 

0.8 

•  2 

0.8 

1 

1.1 

1.1 

1.0 

2.75 

•:  1 

0.9 

1 

1.0 

0.9 

0.8 

2.75 

1 

0.9 

1 

1.1 

1.1 

1.0 

1-2 

0.9 

1 

1.0 

1.0 

0.9 

1-2 

0.9 

1 

1.1 

1.1 

1.0 

•  2 

0.9 

1 

1.0 

1.0 

0.9 

’  2 

0.9 

1 

1.1 

I .  I 

1.0 

Exposure  rate  from  elliptical  cylinder  set 


Exposure  rate  from  circular  cylinder  set 
4  =  Angle  between  receptor-centerline  shortest  distance  line  and  vertical 
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3.2  Concentric  Cylinder  Set  Exposure  Rates 


Each  summation  of  the  basic  equation  presented  above  yields  a  single  exposure  rate 
value,  scaled  to  values  of  unity  for  the  following  parameters:  F,  the  exposure  rate 
constant;  Q',  the  release  rate;  and  u,  the  wind  speed.  To  obtain  corrected  exposure 
rates,  it  is  necessary  to  multiply  the  exposure  rate  values  from  the  curves  by  the 
appropriate  values  of  the  factor  FQ'/u. 

Figure  3B  shows  typical  relative  ground-level  exposure  rates  for  increasing  values  of 
the  ratio  a/a,  which  generally  corresponds  to  increasing  downwind  distance  for  a 
given  set  of  conditions.  This  figure  is  provided  as  an  aid  to  understanding,  and  is  not 
to  be  used  in  calculations.  Figures  4  through  8  provide  the  basic  uncorrected  expo¬ 
sure  rate  data  to  be  used  in  calculations.  The  figures  correspond  to  values  of  a/a 
equal  to  0.3,  0.5,  0.7,  1.1,  and  1.9.  FORTRAN  summations  used  to  obtain  the  data  for 
these  figures  indicate  that  the  maximum  exposure  rate  occurs  for  o/a  -  0.7,  slightly 
less  than  the  maximum  point  indicated  by  differentiation  of  the  dispersion  equation, 
0.707. 

The  ranges  of  parameters  included  in  Figures  4  through  8  are:  5  m  <  L  <  360  m;  E  = 
0.2,  0.5,  and  2.0  MeV,  for  a  -  5,  20,  50,  and  90  m;  E  =  0.5  MeV,  for  a  =  10,  30,  40,  70, 
120,  and  150  m.  All  exposure  rates  are  given  in  mieroroentgens  per  hour.  Interpola¬ 
tion  between  curves  and  figures  may  be  necessary  in  some  eases.  Direct  interpola¬ 
tion  is  often  acceptable;  if  a  more  sophisticated  method  is  needed,  the  Lagrange 
Interpolation  Formula  gives  good  results  (13). 
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3.3  Net  Cylinder  Set  Procedure 


The  basic  CCS  equation  applies  only  to  that  cylinder  set  that  includes  the  receptor  in 
the  plane  of  its  base.  However,  the  correction  for  other  receptor  locations  is 
straightforward  (ref.  10,  p.  369),  and  may  be  accomplished  by  use  of  a  net  cylinder 
set  procedure.  The  steps  in  the  procedure  are  as  follows: 

(a)  Obtain  the  exposure  rate  for  a  set  extending  from  the  plane  of  the  far 
base  of  the  CCS  to  the  parallel  plane  including  the  receptor.  The  corre¬ 
sponding  exposure  rate,  from  this  total  CCS,  is  denoted  -)R(T). 

(b)  Obtain  the  exposure  rate  for  the  CCS  extending  from  th  plane  of  the 
near  base  of  the  CCS  to  the  parallel  plane  including  the  receptor.  The 
corresponding  exposure  rate,  from  this  subtrahend  CCS,  is  denoted  FR(S). 

(c)  Subtract  ER(S)  from  ER(T)  to  obtain  the  exposure  rate  from  the  CCS 
that  is  to  be  evaluated:  i.e.,  the  net  CCS.  This  exposure  rate  is  denoted 
ER(N). 

Values  for  correction  factors,  meteorological  parameters,  and  geometric  parameters 
other  than  cylinder  length  should  be  those  appropriate  to  the  net  CCS,  for  all  three 
sets  used  in  the  net  cylinder  set  procedure.  This  is  an  exact  calculation  procedure, 
not  an  approximation. 

3.4  Application  to  Simple  Flumes 

The  calculation  procedure  for  exposure  rates  applicable  to  simple  plumes  is  presented 
in  Table  4.  The  number  of  cylinder  sets  listed,  including  both  total  and  subtrahend 
sets,  is  sufficient  to  represent  any  simple  plume. 

The  radioactive  decay  factor,  F,  can  be  obtained  by  the  rclutionsh.p  F  =  exp  (-X  x/u), 
with  A  the  radioactive  decay  constant,  s"*;  x,  the  total  path  length  traveled,  m:  and 
u  the  wind  speed,  m/s.  The  eccentricity  correetion  factor,  C,  has  been  discussed 
above.  The  air-ground  interface  scattered  radiation  correction,  A,  will  be  discussed 
in  detail  in  the  discussion  of  irregular  clouds.  The  correction  used  for  simple  plumes 
is  not  as  complicated.  In  general,  a  correction  factor  of  1.2  is  conservative  for  any 
simple  plume,  but  not  overly  conservative. 
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lahlo  4.  Worksheet  Format  lor  Simple  Flume  Calculation 


R  =  X  coordinate  (downwind  distance)  of  receptor 

Multiplied  by  single  scaling  factor  term,  T Q'/u,  applicable  to  exposure  conditions;  +  and  -  indicate  that  exposure  rate 
value  will  be  added  to  or  subtracted  from  sum  of  CCS  exposure  rate  values. 

Total  plume  exposure  rate  is  "sum"  of  14  corrected  exposure  rates.  I.ach  subscript  number  represents  single  value  of  parameter. 


4.0  IRREGULAR  CLOUDS 


The  estimation  of  doses  from  irregular  clouds  involves  assumptions  not  required  for 
evaluating  simple  plumes  (ref.  5,  pp.  89-92).  It  is  assumed  that  the  source  geometry 
of  any  plume  segment  can  be  described  by  use  of  a  plume  trajectory  model  (14). 
Segment  parameters  are  therefore  calculated  as  reflecting  the  plume  segment  path 
history  from  release  point  to  exposure  of  the  receptor.  Figure  9  shows  plume  seg¬ 
ment  path  trajectories  (Figure  9A)  and  calculated  parameters  corresponding  to  those 
trajectories  (Figure  911). 


Figure  10A  represents  the  source-receptor  geometry  of  plume  segments  in  an  irregu¬ 
lar  cloud.  Evaluations  can  be  made  by  visually  inspecting  exposure  rate  curves  for 
(a)  bending  plume  geometry  effect,  (b)  correction  for  anisotropic  diffusion,  and 
(c)  variation  in  effective  stack  height.  The  effect  of  bending  plume  geometry  is  not 
significant  unless  a  plume  changes  direction  sharply  within  a  distance  of  the  receptor 
small  enough  that,  for  that  plume  segment,  the  relationship  L  >  a  holds.  Even  in  this 
case,  the  correction  will  generally  be  on  the  order  of  10%.  A  conservative  dose 
estimate  can  be  obtained  by  using  overlapping  plume  segments  in  the  calculational 
model. 


I  liimc  ID.  I vchnnpics  of it i eiiuhi  1  plume  geometry.  (A)  Determination  <>l'  eyluulci  set  length. 
I.,  ami  teeeptiu-eenteiliiie  distance.  a.  lor  irregular  plume  segments.  (B)  Calculation  ol  ellec- 
t tv o  total  dispetston  coeliivicnl  lot  release  made  under  varying  meteorological  conditions 


Crosswind  anisotropic  diffusion  errors  will  be  relatively  small,  and  the  portion  of  the 
segment  closest  to  the  receptor  may  be  used  to  assign  a  dispersion  coefficient  to  the 
entire  segment  without  excessive  conservatism.  Variable  stack  height  adjustments 
can  be  made  bv  using  different  values  cf  a  for  each  plume  segment,  if  necessary. 


4.2  Calculation  of  Effective  Parameters 


Several  input  parameters  applicable  to  irregular  cloud  volumes  are  calculated  as 
functions  of  the  path  length  traveled  by  a  particular  volume  or  plume  segment.  The 
total  path  length  consists  of  a  series  of  partial  lengths,  or  traverses.  The  calculation 
of  parameters  is  based,  in  some  cases,  on  distance-weighted  averages,  taken  over  the 
sum  of  traverses,  not  the  actual  downwind  distance  of  the  receptor. 

Figure  10B  illustrates  the  method  of  obtaining  dispersion  coefficient  values;  it  in¬ 
volves  incrementing  over  the  sum  of  traverses.  Adjacent  traverses  with  the  same 
meteorological  stability  type  applying  to  them  may  be  treated  as  one  large  traverse, 
for  this  calculation  only.  It  should  be  noted  that  the  minimum  value  of  o  for  a 
particular  traverse  is  equal  to  the  maximum  value  for  the  nearest  upwind  traverse. 
Equations  used  to  calculate  effective  parameter  values  are  presented  below.  In  these 
summations,  the  plume  segment  released  last  will  have  traveled  the  shortest  distance 
and  will  correspond  to  i  =  1,  or  some  other  minimum  value  of  i.  The  earliest  segment 
released  will  correspond  to  i  =  n,  or  some  other  maximum  value  of  i. 

4.3  Equations  Used  to  Calculate  Effective  Parameters 

The  following  equations  are  used  to  calculate  the  effective  parameters  that  require  a 
technique  more  complicated  than  for  evaluation  of  simple  plumes.  All  equations 
apply  to  the  ntfl  traverse  of  a  plume  segment,  and  therefore  include  summed  values 
from  the  entire  plume  segment  path  history  of  that  traverse. 


n 

Total  path  length,  nr.  X,.,  =  X- 

1  i  =  1 


Effective  wind  speed,  m/s:  u  = 


£  <“i  •  xi> 

i  =  1 _ 

n 

V  X. 
i  =  t  ' 


Effective  wind  direction  (*):  w  = 


2J  <wi  ’  xj) 
i  =  1 

n 

I]  X 

i=l  1 


♦Determines  effective  angle  of  segment  centerline.  Use 
any  angular  units,  but  be  consistent. 


Incremented  dispersion  coefficient,  in: 
with:  x.  =  xt  -  x7  and:  a  +  _  = 

1  I  1  Nj. 

Maximum  postulated  segment  length,  m: 


=  S'} 

a  _  _ 

*i  +  i-  si  +  i 


Ln  -  xn  "  xn  -  1  =  un  *  ln 


Definition  of  symbols: 

x  =  distance  from  release  point  to  a  point  downwind,  m 
X  =  length  of  a  traverse,  m 

S  =  Gifford-Pasquill  stability  type,  not  a  parameter 

o  =  increment  or  sum  of  dispersion  coefficient,  m 

u  =  wind  speed,  m/s 

w  =  wind  direction,  angular  units 

x+  =  maximum  downwind  distance  in  a  traverse,  m 

x”  =  minimum  downwind  distance  in  a  traverse,  m 

t  =  length  of  time  applicable  to  a  plume  segment  (release  period  with  constant 
conditions),  s 

L  =  length  of  a  plume  segment,  or  CCS,  m 


4.4  Air-Ground  Scatter  Correction 


The  ratio  of  exposure  rate  at  the  air-ground  interface  to  that  expected  in  an  infinite 
medium  can  be  estimated,  for  irregular  clouds,  by  use  of  an  "image"  source  technique 
(15).  Data  presented  below  have  been  extracted  from  this  published  technique  and 
extrapolated  to  a  wide  range  of  energies  using  other  published  data  (ref.  10,  pp. 
233-245).  The  technique  is  strictly  meant  to  be  applied  only  to  point  sources  but  is 
applicable  to  small  cloud  volumes  located  some  distance  from  the  receptor.  This 
condition  represents  those  cases  in  which  the  technique  is  needed,  since  estimates 
made  for  large  cloud  volumes  located  closer  to  the  receptor  can  be  made  sufficiently 
accurate  by  using  a  simpler  method  applicable  to  simple  plumes,  such  as  assuming  a 
10%  increase  in  the  exposure  rate. 

To  apply  the  technique,  postulate  an  "image"  source  located  at  a  position  in  the 
ground  medium  equivalent  geometrically  to  that  of  the  real  source  (cloud  volume)  in 
the  air  medium.  Assume  that  the  receptor  is  located  on  a  plane  parallel  to  that  of 
the  air-ground  interface.  Refer  to  Figure  11.  Let  hj  and  h2  be  the  minimum  dis¬ 
tances  from  the  receptor  plane  to  the  real  and  image  sources,  respectively.  Let  rj 
and  ro  be  the  distances  from  the  receptor  to  the  sources.  Then  the  angle  V  =  Arcsin 
(h i /r ^ )  is  as  shown  in  Figure  11. 
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The  ratio  of  interfacc-to-infinite  medium  exposure  rates  is  designated  the  air-ground 
scatter  correction  factor,  A.  The  basic  formula  for  calculating  A  is: 


(B 


nr. 


A  =  1  - 


(B  p  )  r; 

/■iTj  2 


1 )  r^  (B  )  r? 

eM(rl"  r2)  •  SA  + - — — —  e^(rr  'V  •  S 


1 


(B  „  )r; 
M  r  j  2 


G 


where  B^r  is  the  infinite  medium  exposure  rate  buildup  factor,  jur  is  the  number  of 
mean  free  paths  between  so  iroc  and  receptor,  and  S.\  and  S(;  are  "air"  and  "ground" 
image  source  strength1',  respectively.  Values  for  these  terms  are  given  in  Table  5. 
Subscripts  1  and  2  in  the  formula  refer  to  distances  between  receptor  and  real  and 
image  sources,  respectively. 


Table  5.  Parameters  and  S()*  of  “Image”  Source  Strength 


Gamma 

Energy 

(MeV) 

Y 

SA 

=  0° 

SG 

H 

*  -C 

15° 

S(i 

SA 

=  30° 

SG 

V 

SA 

-  60° 

s 

f  = 

s\ 

90° 

S(i 

0.20 

0.51 

0.00 

0.29 

0.26 

0.22 

0 . 23 

0.11 

0.27 

0.13 

0.22 

0.50 

0.50 

0.00 

0.27 

0.15 

0.20 

0.20 

0.13 

0 . 20 

0.12 

0 . 1  9 

0.66 

0.49 

0.00 

0.26 

0.14 

0.19 

0.18 

0 . 1  2 

(/.IK 

i).  n 

0.17 

1.00 

0.49 

0.02 

0.23 

0.12 

0.16 

0.14 

0.10 

0.14 

0.  l*i 

(M3 

1.25 

0.49 

0.03 

0.22 

0.12 

0 . 1  5 

0.12 

o.  in 

o .  i  •: 

11  .  1  8 

ii.ii 

2.00 

0.48 

0.04 

0.20 

0,  OH 

0.13 

0.09 

o.ns 

0.10 

o.uii 

•: .  00 

*  S.\  and  S(.  are  terms  for  "image”  source  strength  (ref.  15)  for  air  uni  groins.  . tt- 

buttons,  respectively. 


For  a  receptor  located  at  ground  level,  a  simplification  is  possible,  reducing  the 
calculation  of  A  to  the  form  below.  To  assist  in  calculations  made  with  the  image 
source  technique,  attenuation  parameters  are  given  in  Table  6.  Alternative  refer¬ 
ences  may  be  used  if  desired.  For  a  ground  level  receptor,  use: 

if  >  '  *a  '  s>; 


Table  6.  Attenuation  Parameters  for  Scatter  Correction 

(iamiim  Attenuation  l.xposuiv  Kale  HuiMup  I  actor.  I*  (nr);  r 


fr.nergy 

(MeV) 

(  oof  fie  lent 

10  m 

20  m 

41)  m 

70  m 

10(1  HI 

140  m 

200  m 

0.20 

0.017*8 

1 . 55 

2. 1 1 

3.22 

4.88 

6  .34 

8.78 

1 2 . 08 

o.;>n 

0.0112 

1.20 

i 

1.41 

)  .81 

2.42 

3.03 

3.84 

3.03 

0.68 

Q .0090 

|  1  -  IH 

I  .32 

1  .*54 

2.12 

2.00 

3. 2;> 

4.21 

1.00 

0.0082 

1 .  li 

1.21 

1  .43 

I  .74 

2.11 

2.47 

3.10 

1.2;* 

0.0073 

!  .  09 

1.17 

1.34 

1  .  f»9 

1.87. 

2 .  1  9 

2.89 

2.00 

0.007*7 

i . or* 

1  .  H> 

1  , 20 

1 .30 

1.31 

1.71 

2.02 

.. 

~  -  . . 

- 

.  .  _ 

— 

.  - - . 

loterpol/it  mu  in  ay  Im*  used  for  mt«*rmc«iinte  value's  «.l  r.  or  value**  *  »t  H  obtained  from  'standard  refer¬ 
ences  may  be  "Ubst  itwted. 
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3.0  VARIABLE  RELEASES:  CALCULATION  EXAMPLE 


3. 1  Background  Information 

Postulate  a  release  made  over  a  period  of  50  s  with  a  constant  release  rate  and 
variable  meteorological  conditions,  as  shown  in  Table  7.  For  calculation  purposes, 
the  release  is  arbitrarily  divided  into  ten  plume  segments,  with  each  corresponding  to 
a  parcel  of  gases  released  during  a  5-s  interval. 

A  gamma  level  exposure  rate  estimate  is  taken  at  5  s  after  the  end  of  release  for  a 
location  100  in  from  the  release  point  and  directly  south  (180°)  of  it.  The  elevations 
of  release  point  and  receptor  correspond  to  a  physical  stack  height  of  20  m. 

The  effective  stack  height  for  the  various  stability  types  can  be  calculated  by  use  of 
the  following:  11  =  (20  +  k/u)  m,  with  u  the  wind  speed  in  m/s  and  k  a  constant 
corresponding  to  stability,  with  values  of  8  for  type  B  conditions,  4  for  type  C,  and 

1.2  for  type  D.  The  wind  speed  and  direction  assigned  to  a  plume  segment  are  for  the 
period  between  release  and  dose  estimation. 

5.2  Plume  Segment  Parameters 

Using  the  procedure  described  previously,  calculate  geometric  and  other  parameters 
for  each  of  the  ten  segments.  Geometric  parameters  less  directly  related  to 
meteorological  conditions  are  listed  in  Table  8.  Among  the  parameters  calculated  for 
Tables  7  and  8  are:  (a)  incremented  a,  dispersion  coefficient;  (b)  effective  w,  wind 
direction;  (c)  effective  u,  wind  speed;  (d)  II,  effective  stack  height;  (e)  x,  plume 
segment  downwind  distance;  (f)  v,  receptor-centerline  horizontal  distance;  (g)  a. 
receptor-centerline  distance;  (h)  Lp,  Eg,  and  L^ — the  total,  subtrahend,  and  net 
cylinder  set  lengths,  respectively;  and  (i)  the  ratio  of  o  to  a. 

The  length  of  a  subtrahend  CCS  is  taken  as  the  smallest  difference  between  the 
downwind  distances  of  the  receptor  and  the  CCS.  In  this  example,  the  relationship 
used  to  calculate  x  is:  xj  =  cos  vvj  •  IXj.  I'he  horizontal  receptor-centerline  distance, 
y,  is  taken  at  the  receptor  location,  not  that  corresponding  to  the  CCS.  In  ‘his 
example,  yj  =  sin  Wj  •  100. 

5 . 3  Calculation  Specifications 

Inspection  of  the  data  in  Table  8  and  reference  to  the  CCS  method  permit  the  use  of 
the  following  assumptions:  (a)  If  radioactive  decay  is  insigrificant,  the  contribution 
from  segments  7-10  can  be  ignored  without  underestimating  the  exposure  rate  bv 
more  than  about  5  A,.  This  assumption  requires  the  evaluation  of  both  |xj  -  100|  and 
a.  (b)  For  all  cloud  gamma  dose  equations  of  the  types  estimated  by  the  CCS 
method,  the  exposure  rate  for  all  c,a  <  0.35  can  be  taken  as  equal  to  that  for  a/ a  = 
0.33.  (c)  In  this  example  it  was  assumed  that  E  =  0.5  MeV.  P  =  0.4  R  -  rn-/Ci  -  h, 

and  Q'.  the  release  rate,  is  0.3  mCi/s.  (d)  Correction  factors.  F,  A,  and  C,  arc  taken 
as  a  product  equal  to  1.0.  Based  on  the  above  assumptions,  the  correction  and  scaling 
factor  for  a  plume  segment  is  equal  to  0.12/u,  and  represents  the  product 
I’  •  •  F  •  A  •  C/u.  Table  9  summarizes  the  exposure  rate  calculation  parameters. 
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Stack  Height,  H,  m  24.0  21.3  20.3  21.3  20.2  21.1  20.3  24.0  21.3  24.0 


Table  X.  Geometric  Parameters  of  Plume  Segment 


Plume  Segment  Number 

l 

2 

3 

4 

r' 

6 

7 

8 

8 

HI 

Downwind  Distance  (X  I.  rn 

153 

1-13 

128 

108 

83 

87 

.7 

33 

7T> 

in  ■ 

Horizontal  K-<’  (y(l,  m 

1? 

17 

17 

J  *> 

20 

31 

n 

31 

211 

9  1 

K  -( ‘  UisUim'o  la),  rn 

77 

77 

77 

28 

38 

38 

38 

28 

2*>  | 

Subtrahend.  1. ..  in 

bit 

43 

>8 

8 

7 

,i:i 

48 

t',7 

v:, 

oil  i 

in 

lb 

20 

lb 

77 

1 

2d 

Id 

lb 

Id 

J  oin).  ]  m 

bit 

bH 

48 

24 

48 

'.8 

811 

Hid 

'  a  Ratio,  no  units 

0.48 

11 . 4b 

U.41 

il  .37 

11.31 

II .  ]  “ 

U.  14 

(1.12 

(1.11 

n. nT 

- - - - - _ 

— 

... - 

_ _ - 

_  . 

J 

I  able  c).  Calculation  of  Lxposurc  Rate 


— 

1 

Plume  Seg-  I 
ment  No.  ] 

Ratio, 

Vtt 

CCS  Lengths,  m 

4  4 

— 

a,  m 

I 

i  Sealed  Kxp.  Rate,  fjR/h 

i  EK(T)  LR(S)  EIMN  1 

1  .  .  . 

Corr.* 

Factor 

Corrected 
K.R.,  ^ R/h 

1 

0.49 

i  63 

53 

29 

53 

51 

2 

0.060 

0.12 

2  ; 

0.46 

58 

43 

27 

54 

50 

4 

0.046 

0.19 

3 

0.41 

48 

28 

27 

48 

38 

10 

0.050 

0.50 

4 

0.37 

24 

9 

27 

33 

15 

18 

0.040 

0.72 

5 

0.31 

32 

: 

28 

38 

13 

25 

0.040 

1  .00 

8 

0.30 

48 

33 

38 

i  30 

25 

5 

0.034 

IKK 

_ 

1 

_ 

_ 

2.70 

_ 

•ry'FAC  u 


a . 4  Fxposure  Rate  Estimate 

\ssuming  that  the  total  exposure  rate  is  about  5"<.  greater  than  that  obtained  from 
the  sum  in  Table  9,  the  correct  exposure  rate  can  be  taken  as  about  2.85  £tl\/h.  From 
this,  the  maximum  expected  unrestricted  area  exposure  rate  can  be  taken  as  not 
more  than  25'u  greater  than  the  above  value,  corresponding  to  a  location  at  which 
o  a  0.7.  The  maximum  exposure  rate  is  therefore  estimated  as  about  3.55  /ilt/h. 
and  the  total  exposure  from  the  50 -s  release  is  less  than  0.05  jilt.  Note  that  the  25'b 
assumption  is  valid  only  for  short  downwind  distances;  at  very  great  distances  (c.g.. 
more  than  500  m),  this  assumption  might  lead  to  an  underestimate  of  maximum  expo¬ 
sure  rate. 


SUMMARY 


The  concentric  cylinder  set  (CCS)  cloud  gamma  dose  model  has  been  developed  to 
estimate  external  gamma  doses  from  elevated  radioactive  clouds.  The  primary  appli¬ 
cation  of  this  model  is  for  stack  releases  by  small  radiological  research  facilities, 
including  release  during  changing  meteorological  conditions  or  over  irregular  terrain. 

The  CCS  model  is  applied  with  elementary  calculation  procedures,  using  only  the 
tables  and  working  curves  included  here.  Calculation  procedures  are  provided  for 
dose  estimates  based  on  simple  plumes  and  irregular  clouds.  Sample  calculations  are 
provided  for  evaluating  doses  from  irregular  clouds. 

The  CCS  model  has  been  compared  to  both  published  models,  where  the  latter  could 
be  applied,  and  field  measurements  have  been  made  during  stack  releases  of  argon- 
41,  with  satisfactory  results.  CCS  model  dose  estimates  provide  substantial  informa¬ 
tion  on  environmental  surveillance  for  a  wide  range  of  release  conditions,  with  little 
expenditure  of  time  or  money. 
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